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ABSTRACT 

We have obtained a full suite of Spitzer observations to characterize the debris disk around HR 8799 and to 
explore how its properties are related to the recently discovered set of three massive planets orbiting the star 
We distinguish three components to the debris system: (1) warm dust (T 150 K) orbiting within the innermost 
planet; (2) a broad zone of cold dust (T ~ 45 K) with a sharp inner edge, orbiting just outside the outermost 
planet and presumably sculpted by it; and (3) a dramatic halo of small grains originating in the cold dust 
component. The high level of dynamical activity implied by this halo may arise due to enhanced gravitational 
stirring by the massive planets. The relatively young age of HR 8799 places it in an important early stage of 
development and may provide some help in understanding the interaction of planets and planetary debris, an 
important process in the evolution of our own solar system. 

Subject headings: circumstellar matter - infrared: stars - planetary systems - stars: individual (HR8799) 
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1. INTRODUCTION 

Rapid progress has been made in the past decade in dis- 
covering other planetary systems and in building theories for 
their formation and evolution. More than 300 extrasolar plan- 
ets have been found through measuring stellar radial veloci- 
ties and about 60 planets have been observed to transit their 
stars (Extrasolar Planets Encyclopaedia, http://exoplanet.eu/). 
However, these techniques are optimal in finding systems dra- 
matically different from the solar system, with giant planets in 
orbits very close to their stars. There are two primary avenues 
for probing systems with planets that might have structures 
more analogous to that of the solar system. First, giant plan- 
ets have recently been imaged around Fomalhaut (Kalas et al. 
2008) and HR 8799 (Marois et al. 2008) in orbits with radii 
of tens of AU. Second, the Spitzer Space Telescope {Spitzer, 
Werner et al. 2004) has observed more than 200 debris disks, 
which represent dust production from colliding planetesimals, 
typically in zones analogous to the asteroid and Kuiper belts 
in the solar system. It is expected that these debris disks are 
sculpted by giant planets through gravitational interactions 
(Wyatt 2008 and references therein). Observations and mod- 
els of well-studied examples show a wide variety of structures 
and behavior, possibly related to the individual quirks of their 
(unseen) planetary systems and their evolutionary states (e.g., 
Stapelfeldt et al. 2004; Su et al. 2005; Beichman et al. 2005; 
Song et al. 2005; Rhee et al. 2008; Su et al. 2008; Hillenbrand 
et al. 2008; Backman et al. 2009). 

HR 8799 has been known to have a prominent debris sys- 
tem presumably sculpted by unseen planets (Zuckerman & 
Song 2004). The recent detection of three giant planets or- 
biting this star, along with estimates of their masses and or- 
bits, provides a new opportunity to understand in more detail 
the forces that shape debris disks. This paper reports new 
measurements with Spitzer that fully characterize the debris 
disk around HR 8799. We show that the debris system in- 
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eludes at least three components: warm dust interior to the 
three planets; cold dust exterior to them; and a halo of small 
grains surrounding the cold dust zone. The placement of the 
first two disk components appears to be as expected from the 
orbits of the planets. The dynamics of the three massive plan- 
ets (Fabrycky & MuiTay-Clay 2008; Reidemeister et al. 2009; 
Gozdziewski & Migaszewski 2009) suggests substantial dy- 
namical excitation in the overall planetary system. The ele- 
vated level of dynamical activity in the cold dust zone, which 
is the source of the large halo, is possibly a result of the un- 
stable state of the system. The 20-160 Myr age of the system 
(Mooret al. 2006; Marois et al. 2008) corresponds to a critical 
phase in the history of our own solar system when planet for- 
mation processes were near completion and the solar system's 
dynamical configuration was being established. The HR 8799 
debris system offers an exciting view of this phase of plane- 
tary system formation and evolution. 

We present the new observations in §2, describe the general 
observational results in §3, build detailed models in §4, and 
discuss the implications of this work in §5. 

2. OBSERVATIONS AND DATA REDUCTION 

The observations presented here are from the Spitzer Direc- 
tor's Discretionary Time (DDT) program 530 and Guaranteed 
Time Observations (GTO) program 50175 that utilize all of 
the major observation modes of both the InfraRed Spectro- 
graph (IRS, Houck et al. 2004) and the Multiband Imaging 
Photometer for Spitzer (MIPS, Rieke et al. 2004). The ob- 
servational setting and depth for each of the modes are sum- 
marized in Table 1 . A short (6 s x 1 cycle) integration IRS 
spectrum of HR 8799 was previously published by Chen et al. 
(2006, 2009). Here we present an IRS spectrum that is 4-60 
times deeper along with newly obtained MIPS images. Obser- 
vations at 24 jiva were obtained in standard small-field pho- 
tometry mode, for a total integration of 2 10 s. The observation 
at 70 /im provides a total integration of ^320 s on source. The 
160 /im observation was obtained using the enhanced-mode at 
4 sub-pixel-offset positions, for a total integration of ^150 s. 
The MIPS SED-mode observation provides a low-resolution 
(R= 15-25) spectrum from 55 to 95 /im with 600 s of integra- 
tion on source. 

All of the MIPS data were processed using the Data Anal- 
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Fig. 1 . — MIPS images of HR 8799 with the orientation and beam sizes (FWHMs, black circles) indicated and the surface brightness scale on the side of each 
panel, (a) 24 /xm image before photospheric subtraction, (b) 24 i^m image after subtraction of the model stellar photosphere, (c) and (d) are the 70 and 160 
/im images that are dominated by emission from the disk. The faint 24 irni background source located ~21" from HR 8799 is marked by a black arrow. The 
bright 24 //m source located ~60" away is a field G-type star, BD+20 5278p. 



TABLE 1 
Spitzer Observations of HR 8799 



AOR Key Instrument Module Integration 

28889856 IRS Staring SL2 6 s x 4 cycles 

SLl 14 s X 4 cycles 

LL2 120sx3cycles 

LLl 120 s X 3 cycles 

28889088 MIPS Photometry 24 fim, 5 cluster pos. 3 s X 1 cycle 

28889344 MIPS Photometry 70 fim, default scale 10 s X 4 cycles 

28889600 MIPS Photometry 160 ^m, 4 cluster pos. 10s X 1 cycle 

(enhanced mode) 

25711872 MIPSSED-mode 1' chop 10s x 10 cycles 



ysis Tool (Gordon et al. 2005) for basic reduction (e.g., dark 
subtraction, flat fielding/illumination corrections), with addi- 
tional processing to minimize instrumental artifacts (Engel- 
bracht et al. 2007; Gordon et al. 2007). After correcting these 
artifacts in individual exposures, the final mosaics were com- 
bined with pixels half the size of the physical pixel scale; the 
resulting images are shown in Figure 1 . The calibration fac- 
tors used to transfer the instrumental units to the physical unit 
(mJy) are adopted from the MIPS caUbration papers (Engel- 
bracht et al. 2007; Gordon et al. 2007; Stansberry et al. 2007; 
Lu et al. 2008). 

To ensure the best flat field result at 24 fim, a second flat 
field made from a median stack of the data with all bright 
sources masked out was also applied in addition to the reg- 
ular scan-mirror-dependent flat fielding (for details see En- 
gelbracht et al. 2007). The flux density of the source at 24 
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fxm was estimated using aperture photometry with both small 
and large apertures (the same photometry parameters' as used 
in Su et al. 2006). The large aperture gives a flux density of 
83.0±0.5 mJy, while the small aperture gives a flux density of 
80.8±0.7 mJy (~2.7% lower), implying the source profile is 
more extended than a true point source. Therefore, to capture 
all the flux from the extended component, we conducted pho- 
tometry with a very large (radius of 35") aperture on an im- 
age with the DC offset and all background sources subtracted 
without a sky annulus (i.e., the aperture photometry was rela- 
tive to the median pixel value for the entire frame; therefore, 
an aperture correction of 1.06 was applied, see Engelbracht 
et al. 2007). The final measured flux we adopted is 86.4 mJy 
(before color correction). 

The center of the source at 24 fira was determined using 
a 2-D Gaussian fitting routine, and it coincides with the ex- 
pected stellar position within the pointing error (<1"). This 
is consistent with the fact that the flux contribution in the 24 
/xm band is mostly from the stellar photosphere (58 mJy; see 
§3.1). The steUar photosphere was then subtracted by scaling 
an observed blue point spread function (PSF). After photo- 
spheric subtraction (see Fig. lb), the source has FWHMs of 
5 "78x5. "67 at a position angle (P. A.) of 17°^, slightly broader 
than an observed red PSF (5 "6 1x5. "55 for the ( Lep disk, 
Su et al. 2008). There is a faint (~0.8 mJy) source located 
~21" from HR 8799 at P.A. of -36°, presumably a back- 
ground object. 

The 70 /xm data reduction follows the steps recommended 
by Gordon et al. (2007) and uses time filtering with the source 
region masked out. Several region sizes were tried. A masked 
radius of 55" yields the minimum value for the background 
variation (1-(T7o = 1.14x10"^ mJy arcsec"^ per subpixel). 
The source at 70 /xm (see Fig. Ic) is clearly extended with 
FWHMs of 24"6x24."4 at a P.A. of 77°^ compared to the 
nominal resolution of 18". The source has an azimuthaUy 
symmetric morphology in the image. The outer boundary of 
the source (though masked by the instrumental PSF) can be 
traced out to 42" or 34" in radius at 1- or 3-<tj() levels, respec- 
tively. Aperture photometry was used to estimate the flux den- 
sity for the source since it is extended. An aperture of 42" in 
radius (I-ct boundary) was used with a sky annulus of 44"- 
54". After correcting for lost light based on the same aperture 
settings applied to a theoretical PSF, the final integrated 70 
jim flux density is 545 mJy (before color correction). 

The 160 /um data were taken in the "enhanced AOT" to 
allow them to be time filtered, as was also done for the 70 
fim data (for details see Stansberry et al. 2007). No additional 
160 fim PSF was obtained for the purpose of leak subtraction 
as it has been shown that the ghost image produced by the 160 
/xm filter leakage is less than ~15 times of the photospheric 
flux density at 160 /xm. The disk is ~400 times brighter than 
the photosphere in this channel (see Tab. 2). This is also con- 
firmed by careful inspection of the data where the expected 
ghost image lies. The final 160 /xm mosaic shows some large- 
scale extended cirrus structures surrounding the source (see 
discussion in §5.3). The source is clearly detected at 160 

' smaU aperture: a radius of 6;'23 with sky annulus from 19f'92 to 29;'88 
and an aperture correction of 1.699; large aperture: a radius of 14"94 with 
sky annulus from 29"88 to 42"33 and an aperture correction of 1.142. 

^ The disk brightness in the 24 /xm band is dominated by the bright unre- 
solved component (details see §3.2). Therefore, the position angle given here 
only reflects the angle of the instrumental PSF, not the extended disk. 

' Note that source appears to be azimuthally symmetric, therefore the P.A. 
quoted here does not necessarily reflects the true P.A. of the source. 



TABLE 2 
Observed Flux Densities 



Ac 






Flux Density (mJy) 




l^m 


Total" 


Star* 


Disk MIR Core'' FIR Core'' 


Extended" 


23.68 


86.6 


58 


29 21 <5 


7.5 


71.42 


610 


6.5 


605 • • • 285-387 


319-217 


155.89 


555 


1.4 


554 • • • 554 





"Integrated in a large aperture including the star after color correction 
''From best-fit Kurucz model 

■'Central unresolved source for the inner warm component at 24 /tm 
''Central unresolved source for the outer cold component at 70 fj,m 
"Extended emission for the cold component 



/xm (see Fig. Id), and is consistent with being azimuthaUy 
symmetric above the 5-o-i6o level at the expected star position. 
The source at 160 /im is surrounded by a low level of asym- 
metric background cirrus. The FWHMs are 47"x41" (based 
on an image with a field of view of 88"), slightly broader 
than a nominal observed PSF (~38", Stansberry et al. 2007), 
possibly influenced by the background cirrus. Aperture pho- 
tometry is used to determine the integrated flux. To minimize 
the influence of the background cirrus, we used small aper- 
ture sizes (16" and 24"), sky annuli (64"-128") and aperture 
corrections (4.683 and 2.613 for a red point source of 50 K; 
Stansberry et al. 2007). The resultant integrated 160 /xm flux 
density is 539 mJy (before color-correction). The centroid of 
the disk seen in the MIPS 70 and 160 /xm images coincides 
with the expected position of the star. 

The MIPS SED-mode data were reduced and caUbrated as 
described by Lu et al. (2008) with an extraction aperture of 5 
native pixels (^50") in the spatial direction. Since the source 
size (see below) is smaller than the extraction aperture, the 
slit loss was corrected based on a point source. This proce- 
dure would under-correct for the slightly extended source and 
may explain why there is small offset (~17% at the fiducial 
wavelength of the 70 /im band) between the 70 /im photome- 
try and the MIPS-SED spectrum. The final MIPS SED-mode 
spectrum is smoothed to match the resolution at the long- 
wavelength portion of the spectrum (R=15). 

The basic reduction and extraction of the IRS spectral 
data were provided by the SSC IRS pipeline S18.5. We first 
trimmed a few end points from each module and manually 
removed large outliers by visually inspecting each of the co- 
added, background-subtracted spectra (differencing the two 
nod positions). The spectra were then averaged with sigma 
clipping (S/N=3). All of the final photometry measurements 
(color-corrected, see details in §3.2), MIPS SED-mode data, 
IRS combined spectrum, along with previously published data 
are shown in the spectral energy distribution (SED) (Fig. 2). 

3. ANALYSIS 

3.1. SteUar Properties of HR 8799 

HR 8799 (HD 218396, HIPl 14189), located at 39.4 pc (van 
Leeuwen 2007), is classified as an A5 V (Cowley et al. 1969), 
7 Doradus variable. It has A Bootis type characteristics (Gray 
& Kaye 1999) with low abundances of the heavier elements 
(e.g., [Fe/H] = -0.55), but near-solar abundances of C and O 
(Sadakane 2006). A high-resolution optical spectrum of the 
star indicates T,ff = 7250K, log g = 4.30 (Sadakane 2006). 
The galactic space motion (UVW) of the star resembles those 
of young clusters and associations in the solar neighborhood 
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Fig. 2.— Spectral Energy Distribution (SED) of HR 8799. The best-fit 
Kurucz model (T(,f^=7750 K with log g of 4.5, shown as a grey thin line) 
was determined using ground-based photometry for wavelengths shorter than 
6 /im compiled from the Simbad database (including the 2MASS photometry 
shown as purple filled diamonds, Stromgen photometry shown as red filled 
triangles, Hipparcos BV photometry as blue filled circles, and various John- 
son photometry shown as open squares). The ISO 60 and 90 fim measure- 
ments are from Moor et al. (2006) while the 850 /im datum is from Williams 
& Andrews (2006) and the 1.1 and 1.2 mm observations are from Sylvester 
et al. (1996). The MIPS broad-band measurements are color coiTected. 



with ages of 20-160 Myr (Moor et al. 2006; Marois et al. 
2008). The stellar rotation velocity (vsini) is 49 km s"' (Royer 
et al. 2007), consistent with the star being viewed close to 
pole-on. 

To determine the stellar spectral energy distribution, we 
fit all available optical to near-infrared photometry (Johnson 
UBV, Stromgen uvby photometry, Hipparcos Tycho BV pho- 
tometry, 2MASS JHKs photometry) with the synthetic Ku- 
rucz model (Castelli & Kurucz 2003) based on a good- 
ness of fit test. A value of re//=7500 K with R^=l.4 Rq and 
log g = 4.5 and sub-solar abundances gives the best match, in 
good agreement with the results from spectroscopy (Sadakane 
2006). The resulting stellar luminosity is 5.7 Lq, placing the 
star near the zero-age main sequence consistent with a young 
age. Based on the best-fit Kurucz model, we then estimate the 
stellar photospheric flux densities at 23.6, 71.42, and 155.9 
/xm for the MIPS 24, 70 and 160 fim bands, respectively. 
These values are listed in Table 2. With a stellar mass of 1.5 
Mq, the blowout size (ohi) is ~2 /im assuming a grain density 
of 2.5 g cm-3 around HR 8799. 

3.2. Disk Components 

Figure 3 shows the excess SED after stellar photospheric 
subtraction. It is evident that there are at least two disk com- 
ponents in the system: one with a characteristic temperature 
of ^-^150 K (warm) and the other with a characteristic temper- 
ature of ^^45 K (cold). Similar disk components have been 
suggested by Reidemeister et al. (2009) as well. For black- 
body radiators that are in thermal equilibrium with the stellar 
radiation field, the warm (150 K) and cold (45 K) components 
correspond to radii of -9 AU (0"2) and -95 AU (2."4) from 
the star. Interestingly, HR 8799 d, c and b are at projected 
separations of 24, 38 and 68 AU from the star, respectively 
(Marois et al. 2008). They all lie between the two dust reser- 
voirs. In general, the location of the dust that is at a spe- 
cific thermal equilibrium temperature depends on the grain 
properties (through the absorption efficiency Q„bs, which is 
size dependent). The equilibrium dust temperature is com- 
puted by balancing the absorption and emission energy and 
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Fig. 3.— SED of the HR 8799 Disk (after stellar photospheric subtrac- 
tion). The symbols used are the same as in Figure 2. The short dashed blue 
line shows the blackbody emission at 150 K that matches the excess short- 
ward 30 fj.m well, while the excess longward 30 /im is fit well with modified 
blackbody emission (A^^ '^Ba) at 45 K (long dashed red line). The MIPS 
24 and 70 fim photometry points include the extended component; therefore 
they are slightly higher than the spectra. 



ultimately it depends on cross sections (QahsT^cp-, where a is 
the grain radius); therefore, the temperature of large grains is 
generally lower than that of small grains at the same distance 
from the heating source. We adopt the optical constants for 
astronomical silicates (Laor & Draine 1993) and compute the 
grain properties {Q values for the absorption and scattering 
efficiency for grain radii of 0.1-1000 /im in size) using Mie 
theory with additional modifications for large sizes (a >10 
/im). The resultant thermal equilibrium dust temperatures as 
a function of radius from the star are shown in Figure 4. Due 
to the nature of an optically thin debris system, any dust in a 
continuous spatial distribution between the two distinct tem- 
perature components will have intermediate temperatures, re- 
sulting in a smoother excess SED between 20-55 /im (see 
Fig. 2 in Reidemeister et al. 2009). The SED shape and the 
distinct characteristic dust temperatures strongly suggest that 
very little dust resides between the two components. 

Based on Figure 4, the emitting zone for the warm compo- 
nent is between —5 and ^^15 AU, with an angular diameter 
of < 1". This component, which dominates the emission de- 
tected in the IRS spectrum and the MIPS 24 /im band, should 
be unresolved at 24 /im (inside planet d). A similar calcula- 
tion indicates that the cold component, if it has a size of — 100 
AU in radius, should be at best only barely resolved at 24 /im; 
however, the excess SED (Fig. 3) shows that the cold com- 
ponent would contribute only a small (< 9%) portion of the 
total signal in the MIPS 24 /im band. Thus, the flux well out- 
side the PSF at 24 /im suggests the detection of an extended 
component at this wavelength. 

We estimate the flux of the unresolved point source in the 
24 /xm image (star -i- disk) by scaling a PSF to match the 
peak flux of the source. The scaling suggests that the cen- 
tral unresolved component is —78 mJy at maximum (assum- 
ing the extended component contributes no flux at the posi- 
tion of the star). The stellar photosphere is expected to be 58 
mJy in the 24 /im band, so the central unresolved disk com- 
ponent is —20 mJy before color correction. Due to the large 
temperature differences between the warm and cold compo- 
nents, different color-corrections must be applied to each of 
the components at 24 /im . Applying a color-correction fac- 
tor of 1.055 for a blackbody of 150 K (see the MIPS Data 
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Fig. 4. — Thermal equilibrium dust temperature in the HR 8799 system 
computed based on various grain sizes of astronomical silicates. The two 
horizontal dotted lines mark the dust temperatures (150 and 45 K) derived 
from the excess emission for blackbody radiators. For illustration, the pro- 
jected distances of the three planets are marked as vertical dashed lines. The 
three disk components are marked in different filled patterns (see §3.2 for 
details). 
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Fig. 5. — Radial surface brightness profile of the HR 8799 disk at 70 
/im. Data are shown in open black circles with error bars compared to the 
profile of a PSF (red solid line) scaled to match the peak of the disk flux. An 
extended (outer radius of 46") model surface brightness profile is also shown 
as a dashed blue fine for comparison (For details, see §4.1). The profile at 
radii of 75"-90" is where the second bright Airy ring of the instrumental 
PSF falls, and the match in the profiles between observed point sources and 
the theoretical PSF is poor in this range. 



Handbook, http://ssc.spitzer.caltech.edu/mips/dh/, for the de- 
tailed color corrections for the MIPS bands), the total flux 
density of the central unresolved (warm) disk component is 
mJy, which agrees well with the photosphere-subtracted 
IRS spectrum (22.4±1.4 mJy). Therefore, the extended disk 
component is ~7.5 mJy (after applying a color correction of 
0.894 for a blackbody of 50 K). Adding these measurements, 
the final total flux density for the entire system (star -i- two 
disk components after color corrections for the two different 
temperatures) is 86.6 ±1.7 mJy (assumed 2% error based on 
the nominal calibration uncertainty, Engelbracht et al. 2007). 

The expected stellar photosphere is 6.5 mJy in the 70 
/im band, suggesting the disk is ^603 mJy after a color cor- 
rection factor of 1.12 for a 50 K blackbody. Therefore, the 
total (star -i- disk) flux density at 70 /im is 610±31 mJy (as- 
sumed 5% error based on the nominal calibration uncertainty, 
Gordon et al. 2007). At 160 fim, the stellar photosphere is 1 .4 
mJy; therefore, the total disk flux density is ^^554 mJy after 
a color correction of 1.03 for a blackbody of 50 K, resulting 
in a total (star-i-disk) flux of 555±66 mJy (assumed 12% error 
based on the nominal calibration uncertainty, Stansberry et al. 
2007). 

The 70 fim image is essentially only from the disk since 
the star is <1 % of the total flux in the far-infrared. An ob- 
served radial surface brightness profile of the disk at 70 /im is 
shown in Figure 5. As was indicated by the FWHM mea- 
surements, the disk is clearly resolved. The structure of the 
disk at 70 fim is similar to the one at 24 /im, a bright core 
and an extended disk; however, the bright core at 70 /im can- 
not originate from the warm component seen at 24 /im (be- 
cause the temperature is too cold). This bright core is most 
likely to originate from a dust belt (hereafter the planetesi- 
mal disk) outside planet b, where emission by dust is also 
detected at submillimeter wavelengths (Williams & Andrews 
2006; Sylvester et al. 1996). Consistent with this possibility, 
the emission at 850 /im is mostly confined within 20" (diam- 
eter) based on the JCMT/SCUBA measurements (Williams & 
Andrews 2006). However, if the dust emitting at 70 /im all 
originated from the same location as in the submillimeter, the 
disk could not be resolved by MIPS at 70 /im . In addition, 
Reidemeister et al. (2009) also estimate the outer edge of the 



cold dust component is between 125 and 170 AU (a narrow 
ring) based on simple SED fitting and assumed grain prop- 
erties. Given the nominal resolution of 18" in the MIPS 70 
Hm band, any disk structure that is smaller than 360 AU in 
radius is not resolvable at a distance of 40 pc. The large mea- 
sured FWHMs at 70 /im suggest that the extended emission 
at 70 /im must originate from a component outside the plan- 
etesimal disk. 

Limited by the large beam size at 70 /im, the current data 
cannot differentiate whether there is any gap between the 
planetesimal disk and the extended component. We can only 
estimate the bright core (presumably the planetesimal disk) 
flux by normalizing a PSF to match the peak of the observed 
flux (this procedure provides an upper limit) or to match the 
flux within a 10" radius (to provide a lower limit). After color 
correction for 50 K, the unresolved core flux density is 285- 
387 mJy, suggesting that the extended component is roughly 
equal in integrated brightness to the unresolved core. 

3.3. Disk Inclination 

Based on 10 years of astrometric data on planet b, 
Lafreniere et al. (2009) suggest an orbital inclination of 13°- 
23° off the plane of the sky. In addition, both Fabrycky & 
Murray-Clay (2008) and Reidemeister et al. (2009) suggest 
that the system is unlikely to be exactly face-on from stability 
analyses of the planetary configuration. The resolved image at 
24 /im provides very little constraint because the emission at 
this wavelength is mostly dominated by the central unresolved 
disk. The resolved image at 70 /im appears to be relatively 
azimuthally symmetric. From ^100 MIPS observations of a 
routine calibrator, HD 18071 1 (G9 III, 447 mJy at 70 /im), the 
measured FWHM ratio is 1 .044±0.023 for a point source at 
70 /im. The measured FWHM ratio (1.008) of the HR 8799 
70 /im disk suggests a deviation of ^l.5-a from a face-on 
symmetric disk. From model images (see discussion in §4.2) 
constructed for various inclination angles, we can rule out any 
angles that are larger than ^25° because the resultant images 
give a ratio much larger than the measured value. 

3.4. Summary 
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Fig. 6. — (a) radial profile cut of the HR 8799 disk at 24 fim (black open 
circles). For comparison, similar radial profiles (scaled to match the peak 
of the HR 8799 disk flux) are shown based on the theoretical STinyTim PSF 
(boxcar smoothed, 150 K) and the unresolved disk around HD 38678 (^ Lep). 
The unresolved disk around HD 38678 has high signal-to-noise and shows 
a dynamical range of 2x10* in the radial profile (Su et al. in prep.). An 
extended disk filling the second and third dark Airy rings is evident in the 
profile for HR 8799. (b) radial profile cut of this extended disk at 24 fim after 
subtracting the scaled Lep disk (color-corrected). 



In summary, the resolved disks at 24 and 70 fim indicate 
that the disk around HR 8799 has at least three components: 
a warm belt inside planet d, a planetesimal disk outside of 
planet b with extension of <10", and a surrounding large ex- 
tended disk. The separation into these three components and 
their placement relative to the planets in the system appear 
to be largely independent of assumed grain sizes and optical 
properties. All of the observed flux densities, including those 
for the three disk components, are listed in Table 2. 

4. MODELS 

We now build on the observational results of the preceding 
section through more detailed models of the HR 8799 debris 
disk. Although these specific models are under-constrained 
and hence not unique, they are guided by models of other de- 
bris systems and fit the observations of HR 8799. Therefore, 
they give a plausible picture of the disk behavior In § 4.1, 
we constrain the disk extent by fitting the observed surface 
brightness profiles with assumed powerlaw profiles that are 
independent of grain properties. The main purpose of this 
exercise is to derive spatial information that can be used in 
a realistic thermal model. In § 4.2, we then utilize the disk 
components that are inferred from the observations with the 
adopted grain properties to estimate the physical extent and 
masses required for the dust in each component. 

4. 1 . Disk Surface Brightness Distributions 

Because of our low angular resolution, it is convenient to 
characterize the disk using radially averaged surface bright- 
ness distributions, represented as a mean value at a given ra- 



dius and the standard deviation of the mean. For simplicity, 
we assume that the disk is face-on, so an azimuthally aver- 
aged radial profile can be used for modeling. This approx- 
imation greatly reduces the noise by averaging the intensity 
at a given radius (a similar approach has been used on the 
disks around Vega (Su et al. 2005) and e Eri (Backman et 
al. 2009)). The surface brightness radial profile of the disk 
(after PSF subtraction for the stellar photosphere and mask- 
ing out all nearby sources) at 24 fim is shown in Figure 6a. 
For comparison, we also show the azimuthally averaged ra- 
dial profiles of the smoothed theoretical PSF (generated by 
the STinyTim program, Krist 2006) and the unresolved disk 
around ( Lep"* (Su et al. 2009, in prep.) scaled to match the 
peak flux of the HR 8799 disk. The radial profile of the unre- 
solved ( Lep disk agrees well with the theoretical PSF, while 
it is clear that the second dark and the third bright Airy rings 
of the HR 8799 disk (in the range of -20" to -30" from the 
center) are brighter than for a point source. 

By subtracting the scaled radial profile of the ( Lep disk 
from the observed radial profile of HR 8799, a rough disk sur- 
face brightness distribution for the extended disk component 
is revealed (Fig. 6b). This component is very faint, —1 x 10"^ 
mJy arcsec"^ at a distance of — 1000 AU from the star, similar 
to the surface brightness level seen at 24 fim in the Vega disk. 

As shown in Figures 5 and 6, the true disk structure is 
masked by the instrumental PSFs. One way to understand 
the structure of the disk is to construct a trial surface bright- 
ness distribution, convolve it with the beam, and then com- 
pare with the observed data. In particular, this approach is 
a practical way to constrain the true disk outer radius in the 
low-resolution images. Since our goal is to set constraints on 
the main planetesimal disk, the 24 fim profile in the follow- 
ing fits is for the extended disk only. The simplest model is 
to assume that the surface brightness profile follows a power- 
law, S{r) — r~", between an inner break radius (rbr) and outer 
cut-off radius (ro„,). In an optically thin debris disk where 
the only heating source is the stellar radiation and the disk 
density is only a function of radius in a power-law form of 
— r"'', the radial surface brightness is then proportional to 
r~''B\{Tr) where B\ is the Planck function, and T,- is the radial- 
dependent dust temperature (generally, Tr — r~^-^^ for small 
grains and Tr — r"°^ for large grains). This simple power-law 
formula has been used to fit the well-resolved Vega disk at 
both 24 and 70 fim, and it was found that the power indices 
change from r~^ for the inner part of the disk to r~* for the 
outer part of the disk (Su et al. 2005). An inner break point is 
needed in this simple power-law for steep indices (a —3 or 4); 
otherwise, the profile mimics a point source after convolution 
with the PSF. Because the outermost planet is at a radius of 
—2", we also set the disk surface brightness to zero inside 2", 
i.e., S{r) =0 for r < r/, =2" (h stands for hole). 

We have tried four different cases for the model surface 
brightness distribution: (1) a larger inner hole with a=4, (2) 
a flat distribution between r/, and rbr but with a=4 outward, 
(3) a single powerlaw with a=2, and (4) two powerlaws with 
indices of a=2 and a=3. In all cases, we varied r/„- and ro„, 
of the surface brightness distribution (two free parameters) 
and computed the reduced xl value to determine the best fit. 
The reduced xl, value was computed for data points between 

* Most of the available 24 fim PSF calibrators are stars and thus have very 
blue colors compared to the HR 8799 disk. We use the profile of the Lep 
disk to demonstrate the accuracy of our data reduction and the stability of 
the Spitzer PSF, by matching the theoretical radial profile over a range in 
brightness of over more than 4 orders of magnitude. 
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Fig. 7. — Model and observed surface brightness distributions at 70 fim (red) and 24 fim (blue). Solid and dashed-lines are the model before and after PSF 
convolution, respectively. For details about the model surface brightness parameters see §4.1. In all cases, S(r) = for r <2". (a) is for a large, empty inner hole: 
r,, = r^,. =5"5, S(r) ~ r"^, rout =46"; (b) a flat inner distribution: r;, =2", ri,, =9", S(r) ~ r"", r„,„ =23"; (c) a single power law: r,, =2", rfc^ =10", S(r) ~ r'^, 
Voui =26"; (d) two power laws: r;, =2", Si(r) ~ r"^, rj, =10", S2(r) ~ r"^^ , ro„, =28". The best fit (indicated by lowest xj) is case (a) with Voui =46". 



5"^0", and between 5"-75" at 24 and 70 /xm, respectively. 
Note that the bright hump seen in the 70 /im profile between 
75" and 90" is where the second bright Airy ring is located 
(see Fig. 5). The match between an observed PSF derived 
from the calibration stars and theoretical PSF is poor in this 
range; therefore, we disregard the 70 /im data points outside 
75" in the modeling. The results are shown in Figure 7. 

Due to Spitzers large beam sizes in the far-infrared, the as- 
sumed surface brightness behavior and its number of free pa- 
rameters needed to describe the model, there is no unique fit 
for the observed profiles. The best-fit outer radii range from 
23" to 46", consistent with the previous finding of emission 
from outside of the submillimeter planetesimal ring (<10"). 
From the model surface brightness profiles, we can also de- 
rive the surface brightness ratio {870/824) in the extended disk. 
The ratios range from ^170 to ^-^250, which can be used 
to set constraints on the grain properties since the density- 
dependent term cancels out in the ratio. Figure 8 shows the 
expected ratio for a size distribution of a~^-^ with a maximum 
size cutoff of '^lO fim and various minimum grain size cutoffs 
(amin,cutoff) using astronomical silicate grains. Based on the 
range of derived 870/824 ratios (represented by the horizontal 
dashed line in figure 8), the acceptable range of amm,cutoff is 
<2 /im {=ai,i). Outside the main planetesimal disk (r >10"), 
only grains with minimum sizes < the radiation blowout size 
have flux ratios in the observed range. This supports the idea 
that the disk outside of the main planetesimal disk consists of 
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small grains being pushed outward by radiation pressure, i.e., 
the HR 8799 disk has an extended halo similar to that seen 
around the Vega disk (Su et al. 2005). 

4.2. 8pecific Disk Models 
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Fig. 9. — Model SED (top panel) and surface brightness profiles (bottom panel) compared to the data. Two sets of models are shown having the same 
inner warm disk component warm =6 AU, Roui .warm='\-5 AU, flm/n=l-5 ^m, amax=^-5 /^m with total Md of 1.1 x 10~^M^) and (a) a narrow planetesimal disk 
model {Rm,pb=90 AU, Roiir,pb=^50 AU, a,nin='\.0 fim and a,„a.x=lOOO fim with M(/=5.5x lO^^M^) plus a halo (RinMio=^50 AU, Rouij,aio=^000 AU, £/«„„=l 
fim and a,nax=^ A^ni with Md=\.2x\0~^M(^), or (b) a wide planetesimal disk model (7?,>,^fo=90 AU, Rf,i,r pij=300 AU, a,ni„=\0 //m and a,i,ax=^000 fim with 
Mrf=1.2x 10"^M©) plus a halo (^,„./,«/<,=300 AU, /?„„,./„,/^=10()0 AU, <3„„„=1 fim and (3,„«.v=10 //m with Mrf=1.9x IQ-^M©). The narrow ring planetesimal disk 
model gives a satisfactory fit to the SED, but fails to match the disk radial profiles in contrast to the fit of a wide planetesimal disk shown in (b). The model 
submillimeter flux appears to be too high compared to the observed data. However, because of the source extension, the "observed" model flux actually agrees 
with the measurements that were calibrated assuming an unresolved source. For details see §5.1.2. 



From the preceding section, we can expand on the descrip- 
tion of the three disk components summarized at the end of 
§3.2: (1) the warm beh inside the innermost planet d has a 
temperature of '--^150 K with parameters dubbed as "warm"; 
(2) the main planetesimal disk outside the outermost planet b 
has a temperature of ~45 K with parameters dubbed as "pZ?" 
(parent bodies); and (3) the outermost component is small 
grains strongly affected by photon pressure with parameters 
dubbed as "halo". The basic equations for computing the 
SED and surface brightness of an optic ally- thin debris disk 
can be found in Backman & Paresce (1993), Wolf & Hillen- 
brand (2003), and Su et al. (2005). To model these structures, 
we use grain properties for astronomical silicates with a den- 



sity of 2.5 g cm"-' and a n{a) ^ a~^'^ size distribution appro- 
priate for a steady-state collisional cascade (Dohnanyi 1969). 
Using the SED of the excess and the resolved images as con- 
straints, we fit the minimum and maximum grain sizes (flm,„ 
and Umax) in the size distribution and the inner and outer radii 
of each of the disk components (7?,„ and Rom) by assuming 
a constant surface density (S(r) ^ r") disk for the strongly 
bound disk and S(r) ~ r~' for the weakly-bound or unbound 
disk (i.e., the grains affected by photon pressure). 

From its characteristic temperature (^^150 K) and Figure 4, 
the inner radius of the warm component can be as close as 
~6 AU if the grains have a size of '-^lO /im. However, the 
spectral shape of the 10 to 20 fim excess emission is better 
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fit with grains of ^fim sizes. In addition, the spectral shape 
suggests that sub-micron grains are absent, since their emis- 
sion spectrum would result in higher contrast features at 10 
and 20 fjm. We fit the inner warm component as a bound disk 
with parameters used in Table 3. The dust mass only accounts 
for the small grains that dominate the mid-IR emission. We 
lack the constraints, such as the flux of this warm component 
at longer wavelengths, that would be required to estimate the 
mass of the population of large grains and parent bodies asso- 
ciated with the warm emission. 

It is less straightforward to fit the other two disk compo- 
nents, due to the lack of strong constraints on the size of the 
planetesimal disk. We can, however, use the following infor- 
mation: (1) we know roughly how much flux comes from the 
unresolved component (presumably the planetesimal disk) at 
70 /im; (2) the minimum inner radius of the planetesimal disk 
is ^90 AU, from the characteristic temperature of the cold 
component in the SED; and (3) large grains (at least a few 
hundred microns) are required due to the coUisional cascade 
nature of debris disks. The model parameters not only have 
to provide satisfactory fits to the global SED, but the model 
images after convolution with the beam sizes have to match 
the observed images as well. Our modeling strategy is similar 
to the one applied to the 7 Oph disk (Su et al. 2008); we first 
tried many combinations of parameters that can provide good 
fits to the SED, then refined the parameters by matching the 
observed radial profiles. We started with constraining the pa- 
rameters for the inner component using data shortward of 30 
fjxa since they are most sensitive to this component, then used 
the constraints Usted above for the parameters in the other two 
disk components. 

For the planetesimal disk, we fixed and Ri„,ph at 1000 
/zm and 90 AU, respectively. Additionally, we constrained 
amin to be greater than a^/ (~2 ^m). For the extended halo 
disk, RoutMio was fixed at 1000 AU; Rm^haio and Rout,pb were 
required to be equal (there is no gap between the planetes- 
imal and halo disks). Initially, we fit the planetesimal disk 
with a narrow ring (A/? < R) so the actual density distribu- 
tion has less effect in the output SED. A difficulty in this nar- 
row ring model is that grains smaller than au have relatively 
high thermal equilibrium temperatures at the location where 
the extended disk starts, and the resultant emission from the 
inner edge of the halo component is too strong in the 25-35 
jiva range. The only way to make the narrow ring model 
work is to have far less contribution from the extended disk 
as shown in the SED of Figure 9a (the top panel). However, 
the resulting model yields insufficient flux from the extended 
component, and therefore, also shows a disk at 70 //m that is 
too small, as shown in the model radial profile (Fig. 9a). 

To satisfy all the constraints listed above simultaneously, 
the outer radius of the planetesimal disk has to extend to 
^^300 AU. Figure 9b shows one of the best-fit models in terms 
of SED and radial profiles. In this model, the planetesimal 
disk extends from Ri„,pi,=9Q AU to Rom ,pi,=300 AU and con- 
sists of grains of amm=^0 /um to aTO„=1000 /um with a total 
dust mass (M^) of 1.2x 1O"^M0. The halo disk extends from 
Rin,haio=^00 AU to ^?oM/,/!a/o=1000 AU and consists of grains 
of flm,„=l /xm to 0^01=10 um with a total Mrf=1.9x lO'^M^. 
Note that the reduced xt values in the SED fitting are similar 
between the narrow-ring and broad-disk models; however, the 
model flux ratio at 70 /xm between the extended component 
(halo) and the unresolved core (planetesimal disk) is too low 
compared to the observed ratio (~1). This mis-match is best 



TABLE 3 

Parameters in the Preferred Model 



Parameters Inner Warm Disk 


Planetesimal Disk 


Halo 


(warm) 


(pb) 


(halo) 



S(r) 




^^0 




R,„ (AU) 


6 


90 


300 


Rou, (AU) 


15 


300 


1000 


Omin (Aim) 


1.5 


10 


1 


amax (/im) 


4.5 


1000 


10 


Ma (Me) 


1.1x10-" 


1.2x10-' 


1.9x10-2 


fd = Lir/L^. 


2.2x10-5 


7.5x10-5 


1.0x10-* 



shown in the model surface brightness profile at 70 /im (bot- 
tom panel of Fig. 9a), where the narrow-ring model does not 
fit the 70 /xm profile at all. 

The parameters used and derived in our thermal model for 
the disk components are listed in Table 3. Because of the am- 
biguous nature of SED modeling with multiple components 
and low-resolution images, these model parameters are not 
unique but were built from the least number of components 
with simple assumptions that produced a good match to all 
the available data. One should not take the model grain sizes 
too literally since they depend greatly on the grain properties 
used. However, the total dust mass in each of the components 
remains similar among aU of the acceptable models. 

5. DISCUSSION 

5.1. Debris Disk Structure 

The architecture of the HR 8799 planetary system is com- 
plex: a warm inner asteroid-belt analog located at a radius of 
6-15 AU; a cold Kuiper-belt-analog planetesimal disk located 
from ~90 AU up to ~300 AU; three massive planets orbiting 
between the two disk components; and a prominent halo of 
small grains extending up to ~1000 AU. Similar components 
(except for imaged planets) have been found in a number of 
other debris disk systems, indicating some underlying order 
within the diversity in debris disk structures. 

5.1.1. Inner Warm Component 

Assuming the grains in the inner belt are astronomical sili- 
cates, we can rule out grains of sub-micron size, due to the 
absence of strong silicate features at 10 and 20 /xm. The 
shape of the IRS spectrum strongly favors micron sizes for the 
grains that dominate the emission. If we further assume that 
all the grains in this component are bound (am,„ = ai,i), then 
the outer radius of this warm component can be as smaU as 
~ 10 AU based on temperature arguments (see the 2 /xm curve 
in Fig. 4). We also find that the maximum inner radius of this 
component is around ^6 AU, again from temperature argu- 
ments. There are few, if any, grains inside this radius. An 
upper limit on the dust mass inside this inner warm compo- 
nent can be estimated based on the observed IRS excess spec- 
trum. A flux density of 1 mJy at 10 /xm corresponds to a 
dust mass of 2.3x lO'^M© assuming astronomical silicates of 
fl=10 /im at 3 AU (i.e., T4 ^200 K) and a grain density of 2.5 
g cm"-'. This inner hole can be maintained by ice sublimation 
if the grains are icy (Jura et al. 1998); a process that will oc- 
cur at a dust temperature of ~150 K. Alternatively, the lack of 
grains inside ^6 AU suggests either that there is another (un- 
seen) interior planet maintaining a dust-free zone, or that the 
growth of the protoplanets in the terrestrial zone (3-20 AU for 
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an A-type star) at an early stage initiated a collisional cascade 
that moved outward and created this inner hole (Kenyon & 
Bromley 2004). Since we lack constraints such as a measure- 
ment of the long wavelength emission from this inner com- 
ponent at higher spatial resolution, we cannot set limits on 
Umax- Because the dust in this component is warm, very lit- 
tle mass is required, only about 1.1 x lO'^'M® for dust with an 
infrared fractional luminosity (fd) of 2.2 x 10"^. These values 
are rough lower limits because it is difficult to place tighter 
constraints without knowing additional parameters such as 

^max- 

We also fit the inner warm component with amorphous car- 
bon grains (density of 1.85 g/cm^, Zubko et al. 1996). As 
expected from the discussion in § 4, the model results for 
the placement of the disk are not very sensitive to the as- 
sumed grain composition warm ~ 6 AU and Rout,warm ~15 
AU), and we can also rule out a significant amount of sub- 
micron carbonaceous grains in this inner component because 
the equilibrium temperature of 150 K places them right at the 
location of planet c. A similar result of low sensitivity to grain 
properties is also found in Su et al. (2005), who analyzed the 
behavior of the Vega system for grains of silicate, carbona- 
ceous, and mixed composition. 

A similar inner warm component (T^ ~150 K) that dom- 
inates the disk emission at the IRS and MIPS 24 fim wave- 
lengths is seen relatively commonly among resolved de- 
bris systems. For example, HR 4796 A, Fomalhault, and e 
Eri disks all have such a component (Wahhaj et al. 2005; 
Stapelfeldt et al. 2004; Backman et al. 2009). In fact, all 
of their excess SEDs look very similar: a steep inflection 
between the IRS spectrum and MIPS-SED mode data, indi- 
cating a relatively sharp inward-facing edge to the structure 
dominating the far-infrared emission. This feature may arise 
from the action of massive planets, as appears to be the case 
for Fomalhaut (Kalas et al. 2008; QuiUen 2006; Chiang et al. 
2009). HD 32297 and r? Tel are other systems that show sim- 
ilar characteristics revealed by recent ground-based imaging 
(Moerchen et al. 2007; Fitzgerald et al. 2007; Smith et al. 
2009). In addition, a recent IRS spectroscopic study of un- 
resolved debris disks shows that warm components are very 
common (~50%) among A-type excess stars (Morales et al. 
2009). In the case of HR 8799, there are strong indications 
that this warm component is physically separated from the 
cold, outer component and that there is little dust in the inter- 
vening region. 

5.1.2. Outer Cold Planetesimal Disk 

Our best-fit model appears to indicate too much flux in the 
submiUimeter (Fig. 9) for the cold planetesimal disk. How- 
ever, the model 850 /xm disk has a FWHM of 17."8 after being 
convolved with the nominal 850 fjm beam size of 15". This 
slight extension would not be resolved by the 850 obser- 
vation because the data were taken in photometry mode and 
no mapping was done. The flux density reported in Williams 
& Andrews (2006) was calibrated assuming a point source. 
Therefore, the correct way to compare the model and ob- 
served fluxes is to integrate the total flux in a 15" diameter 
aperture on the model image. The 850 /im model flux within a 
diameter of 15" is 9.33 mJy, consistent with the SCUBA mea- 
surement of 10.3±1.8 mJy. Similarly, the 1.2 mm flux is also 
under-estimated for the source extension since it was obtained 
using the IRAM 30 m single dish telescope (FWHM=9"5). 

The dust mass in the planetesimal disk that is derived from 
our three-component model is 0.12 Mq, consistent with the 



mass derived from the 850 /im observation (0. 1 M®, Williams 
& Andrews 2006). This is only twice as large as the mass re- 
quired by the narrow ring model, even though the disk is much 
larger. The dust mass in this cold component is similar to val- 
ues derived for disks in the 12 Myr-old (3 Pic moving group 
(/3 Pic and HD 15115, Holland et al. 1998; Williams & An- 
drews 2006). The amount of the excess emission in the HR 
8799 disk (/?24= 1 -5 and R-io=9A, compared to the stellar photo- 
sphere) is comparable with other debris disks around A-stars 
of similar age (20-160 Myr) (Su et al. 2006), at least within 
statistical expectations. The infrared fractional luminosity for 
this component is 7.5 x 10"^, suggesting a collision-dominated 
disk. 

Since the cold component is the birth place for the second- 
generation debris from coUisional cascades, the presence of 
large grains is required. The inner edge of the disk, although 
not directly resolved by imaging, has a minimum radius of 
~90 AU from temperature arguments. We lack high-spatial- 
resolution submiUimeter data to put strong constraints on the 
outer edge of the planetesimal disk. The outer radius of ^^300 
AU in our model is based on the fewest and most straightfor- 
ward assumptions. 

5.1.3. Halo 

Finally, we consider the halo component of the system. The 
general appearance and extent of this component are similar 
to the halo around the Vega system. In that case, the Spitzer 
observations provided sufficient angular resolution to fit the 
surface brightness profile and show that it is consistent with 
expectations for grains being blown out by radiation pressure. 
While the details of this process are complex, it is generally 
parametrized by (3, the ratio of radiation pressure force to 
gravitational force. When collisions produce grains with (3 
> 1, they will be blown out of the system quickly. When the 
collisional cascade produces grains over a narrow range of /? 
below unity, they will tend to go into elliptical orbits. The 
behavior of these barely-bound grains has been suggested to 
explain the broken power-laws seen in the scattered-light disk 
surface brightness profile around AU Mic (Strubbe & Chi- 
ang 2006) and (3 Pic (Augereau et al. 2001). However, these 
highly-elliptical orbits will continue to cross the parent body 
zone and eventually it is likely that these grains will collide 
again, be broken into smaller ones with (3 > \ and be blown 
out of the system. Thus, the time scale for mass loss from the 
system depends on the details of the coUisional fragmentation 
of the grains. 

However, qualitatively, the halo around HR 8799 resembles 
the behavior of Vega closely and differs substantially from the 
observations of most other A-stars, indicating that both sys- 
tems are currently undergoing an enhanced level of dynami- 
cal stirring in their outer planetesimal disks that leads to an 
elevated production rate for small grains with (3 near unity. 
These grains will eventually be ground down to sizes that par- 
ticipate in a general outflow from the system. We therefore 
model the halo by making assumptions similar to those used 
for the Vega outflow (Su et al. 2005). 

To fit the Vega outflow, grains extending up to about five 
times the nominal blowout size were required. This apparent 
discrepancy probably arises because, for simplicity, the mod- 
els are based on solid, spherical grains, whereas real grains 
are Ukely to be non-spherical and complex in structure with 
a substantial volume of voids (e.g., Dominik et al. 2007 and 
references therein). The influence of photon pressure will be 
much larger for such grains than for solid spherical ones. We 
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use a„,ax =10 /im in the extended halo component, about five 
times larger than the blowout size for HR 8799, and other pa- 
rameters of the model are similar to those used for Vega (Su 
et al. 2005). The estimated dust mass in the halo is relatively 
large, 1.9x lO'^M^ that is ^^6 times greater than the mass in 
the Vega outflow. The terminal radial velocity is 2-A km s"' 
for grains having (3^1 and that are being produced in the 
planetesimal disk, 100-300 AU from a 1.5 Mq star (Moro- 
Martm & Malhotra 2005). For a disk radius of 1000 AU, the 
residence time for such particles is ^ 900-2400 yr. For a 
static collisionally dominated disk (Wyatt et al. 2007), the ex- 
pected ratio of mass in blown out grains to that of grains in the 
bound disk is -1.2% in the HR 8799 system (i.e., a 1.5 Mq 
star with a fractional disk luminosity of 8 x 10"^ and r '--^ 100 
AU). The observed ratio for HR 8799 is 15 times higher, qual- 
itatively similar to the Vega system, where the mass in blown 
out grains is about 30 times the predicted static level. Both 
systems have much more prominent halos than systems such 
as Fomalhaut (Stapelfeldt et al. 2004) and (3 Leo (Stock et al. 
in prep.). This comparison indicates that there is an elevated 
level of dynamical activity in the bound disks of Vega and HR 
8799 that is enhancing the production of small grains. 

Although we have proceeded by analogy with Vega, future 
work should address the role of weakly bound grains in more 
detail. The conclusion that the halo indicates more dynamical 
activity than in A stars without such halos should be robust, 
but more complete modeling should draw out other aspects of 
the HR 8799 system. 

5.2. Interaction of Planets and Debris 

Based on the "resonance overlap" condition in the planar 
circular restricted three body problem, the width of the chaotic 
region, Aa, in the vicinity of a planet with mass of Mp and 
orbital semi-major of ap orbiting a star with mass of is 
given by Malhotra (1998), 

Afl~ 1.4fl^(M;,/M,)<2/7>. (1) 

Adopting the nominal parameters from Marois et al. (2008), 
stellar mass 1.5 Mq, planet masses 7 Mj and 10 My, and semi- 
major axes 24 or 68 AU for planets d and b, respectively, the 
width of the unstable zone near each planet is is —8 and —23 
AU, respectively. Therefore, the locations of the outer edge of 
the inner dust belt (~15 AU) and the inner edge of the outer 
belt (—90 AU) are both consistent with gravitational sculpting 
by the innermost known planet d and the outermost known 
planet b, respectively. The nearly dust-free zone interior to 
—6 AU may indicate the presence of additional unseen inner 
planets. If so, then this leads to an estimate of the timescale of 
inner planet formation of <20-160 Myr given the age of HR 
8799. This is in good agreement with theoretical calculations 
for terrestrial planet accretion (Wetherill 1992) as well as iso- 
topic constraints on the age of the Earth (Podosek & Ozima 
2000). 

The existence of the halo indicates that the outer planetes- 
imal disk is heavily stirred due either to: (1) the planets 
not being in a stable configuration (Fabrycky & Murray-Clay 
2008; Reidemeister et al. 2009; Gozdziewski & Migaszewski 
2009); (2) planet migration and orbital resonance crossings 
that are causing extreme excitation in the planetesimal disk, 
as in models for the Kuiper Belt (Malhotra 1995) and for the 
Late Heavy Bombardment (Gomes et al. 2005); or (3) pro- 
cesses in the planetesimal disk itself, such as ongoing forma- 
tion of (unseen) ice giants or planetary embryos that stir the 
disk (Kenyon & Bromley 2008). The first two possibilities 




Fig. 10. — The CO J=3-2 contours (green) overlaid on the 160 data 
(pseudo-color with blue contours). The field of view is 7' x5.9' with N up 
and E toward the left. The contours of the 160 /jm data range from 1-, 2-, 
3-, 5- and S-cr detection level, while the green CO contours show 5 linear 
levels of integrated intensity from 0.94 (2-cr) to 2.7 (7-(t) K km s"' . The 
stellar position is marked as a star-shaped symbol and the blue aiTow shows 
the proper motion direction of HR 8799 with a length indicating 40" . 



relate the dynamical activity in the planetesimal disk to is- 
sues associated with the stability of the configuration of the 
three known planets, while the third hypothesis would be in- 
dependent of the planets. Although no specific conclusions 
can be drawn, the outcomes of any of these scenarios are a 
heavily stirred planetesimal disk that produces an era of en- 
hanced collisional activity. The resulting collisional cascades 
and avalanches will feed an outflow of small grains (Grig- 
orieva et al. 2007). 

5.3. The Background Cloud near HR 8799 and the X-Bootis 
Phenomenon 

The cirrus seen at 160 fim resembles some of the large-scale 
structure in the CO (3-2) cloud detected by Williams & An- 
drews (2006). Figure 10 shows a new, 7' x7' CO (3-2) map at 
15" resolution around HR 8799, taken with the HARP hetero- 
dyne array on the JCMT, overlaid on the Spitzer 160 fim map. 
The CO emission is integrated in the local standard of rest 
(LSR) velocity from -4 to -6 km s"'. The radial velocity of 
the star, -12.4±0.5 km s~', from Moor et al. 2006 is in the 
heliocentric velocity frame, however, and converts to -4.6 km 
s"' LSR. Thus, contrary to the conclusion in Williams & An- 
drews (2006), we cannot say that the cloud is not associated 
with HR 8799 and cannot place any limits on the gas content 
of the disk from these observations. This is also further con- 
firmed by the new IRAM 30-m observation that found strong 
CO (2-1) line emission due to the cloud at the heliocentric 
radial velocity of HR 8799 (Hily-Blant, Kastner, Forveille & 
Zuckerman, in prep.) 

This new JCMT CO (3-2) map is larger and of higher qual- 
ity than the map presented in Williams & Andrews (2006) and 
shows the extended, low level molecular emission in more de- 
tail. It appears to be associated with the extended arc of high 
latitude clouds, MBM 53-55, detected in CO (1-0) by Mag- 
nani et al. (2000). The distance of these objects and their rela- 
tion to one another is not well known. If the association with 
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the star is physical then, at 40 pc, this would be the closest 
known molecular cloud. 

Whether physically associated or not, it is hkely that the 
low flux levels (l-3o-) of the 160 fjm emission come from the 
molecular cloud. The asymmetry in the outer part of the ob- 
served source at this wavelength may also be due to cloud 
contamination. Nevertheless, the detected source (with the 
peak flux > 8-(t) is coincident with the expected stellar posi- 
tion, which is offset (~20") from the central emission of the 
nearby clump in the cloud. The cloud therefore contributes 
only a very small amount of flux that cannot be completely 
removed from our aperture photometry. At 70 fim, there is 
no resemblance to the large-scale background structure in the 
data nor an asymmetry seen in the source morphology; there- 
fore, we can rule out any significant contribution from this 
extended background cloud to the 70 /zm disk emission. 

The arrow in Figure 10 shows the direction of proper mo- 
tion of HR 8799 (SRA = 107.93 mas yr"' and SDec = - 
49.63 mas yr~'; van Leeuwen 2007). Given the coincidence 
in space and velocity, it is plausible that accretion from the 
cloud might explain the abundance anomaUes leading to the A 
Bootis designation (Kamp & Paunzen 2002). Using equation 
(16) from Martinez-Galarza et al. (2009) for Bondi-Hoyle ac- 
cretion implies a mass accretion rate of ~ 3 x 10~^'^MQyr~^ 
with an assumed density of 100 cm~^, just above the thresh- 
old to produce A Bootis behavior (Turcotte & Charbonneau 
1993). However, this calculation assumes a fluid model for 
the ISM, which is not strictly applicable at the expected den- 
sity (Martmez-Galarza et al. 2009). Therefore, the accretion 
rates may be significantly lower than indicated (Alcock & II- 
larionov 1980). The accretion hypothesis rests on the untested 
assumption that the rate can be increased through such means 
as small scale high-density structures in the ISM and/or in- 
creases in density in the wake of the star. 

6. CONCLUSIONS 

We have obtained deep infrared images and spectra of the 
debris system around HR 8799. To our surprise at the distance 
of ~40 pc, the HR 8799 disk is clearly resolved at both 24 and 
70 /xm; the outer boundary of the disk can be traced further 
than 1000 AU. From general, model-independent considera- 
tions, the star appears to have one zone of warm dust within 
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the orbit of its innermost known planet, another broad zone of 
cold dust outside the orbit of its outermost known planet, and 
an extended halo of small grains from this outer zone. 

We have constructed detailed models of the system, guided 
by models of other debris systems as well as the observations 
of HR 8799. Although the models are not unique, because we 
draw on experience in modeling many debris systems, they 
should give a reasonably accurate picture of conditions in the 
HR 8799 system. The models indicate that the inner zone lies 
between 6 and 15 AU with a total dust mass of 1.1 x lO'^M^ 
and the outer zone is between 90 and 300 AU with a dust 
mass of 1.2x1O~'M0. There is little dust inside the inner 
zone (r <6 AU) and in between the inner and outer cold zones 
where the three planets reside. The halo extends to at least 
1000 AU and has a dust mass of 1 .9 x 10~^M^, about 15 times 
the level expected for a static, quiescent debris system. The 
implied high level of dynamical activity in the HR 8799 debris 
system may be related to perturbations imposed by its system 
of three massive planets or may be signaling ongoing planet 
formation in the outer parts of this system. 

HR 8799 is at an age similar to that of the solar system 
when the terrestrial planets formed, along the way to settling 
into the final configuration. The planetary system of HR 8799 
(both planets and debris disk) provide an intriguing snapshot 
of processes occurring at this stage and may help us under- 
stand the formation and evolution of our own planetary sys- 
tem. 
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